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Abstract 

We continue our investigation of massive gravity in the massless limit of vanishing gravi- 
ton mass. From gauge invariance we derive the most general coupling between scalar matter 
and gravity. We get further couplings beside the standard coupling to the energy-momentum 
tensor. On the classical level this leads to a further modification of general relativity. 
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1 Introduction 



In a previous paper [T] the massless limit of massive spin-2 quantum gauge theory (called massive 
gravity for short) has been studied, and the corresponding classical theory has been taken as an 
alternative to general relativity. The reason why the limit of vanishing graviton mass m — > 
leads to a modification of general relativity is the vector graviton field which is essential in 
the massive theory and does not decouple from the symmetric tensor field h^^ in the massless 
limit. However, the theory studied in [l] is not yet complete because the coupling to normal 
matter was described in the standard way by means of the energy-momentum tensor of ordinary 
matter. Then there is no direct coupling between the f-field and ordinary matter. But it is our 
ultimate aim to derive all couplings from gauge invariance. For the scalar matter couplings this 
is done in this paper. We will find that further couplings between the scalar field $ and the 
vector graviton field with arbitrary coupling constant are possible in massive gravity. 

As in [1] we then consider the massless limit m ^ 0. In order to get a non-trivial limit we 
have to choose the free coupling constants proportional to the graviton mass ra. This is not 
unusual because coupling terms with factors m appear also at other places in massive gravity. 
The surviving couplings then lead to additional terms in the classical Lagrangean. This is a 
further modification of general relativity. 

The paper is organized as follows. In the next section we derive the most general gauge 
invariant trilinear coupling between scalar matter and massive gravity. We apply the descent 
method which was already used for the construction of pure massive gravity [7j. To get unique- 
ness of the result the cohomological methods developed in [3] and [5] have to be employed. We 
find five possible couplings where three contain the vector graviton field v^^. 

In Sect. 3 we study second order gauge invariance which gives further restrictions on the 
coupling. Only three coupling terms survive: one is the well-known coupling to the energy- 
momentum tensor of the scalar field <I>, the second is the self-coupling and there is one 
coupling to the w-field. However, the necessary finite renormalizations generate new quartic 
couplings. 

In the last section we investigate the new couplings in the limit of vanishing graviton mass 
m ^ 0. Since two of the quartic couplings contain m in the denominator, a non-trivial limit 
only exists if the (free) coupling constants are proportional to m. The new coupling terms then 
give rise to the modification of general relativity mentioned above. 

2 Gauge invariant couplings to scalar matter 

In [6] we have analyzed the interaction of massless gravity with massive Yang-Mills fields and 
with scalar fields. The coupling of the free quantum fields can be obtained with the cohomology 
methods developed in [4J and [5j. The case of massive gravity can be analyzed with the same 
methods. We work in the same setting as in [l], in particular we have the gauge structure on 
the free asymptotic fields defined by the gauge charge operator Q and the corresponding gauge 
variation dq. First we give the expression of the gauge invariant variables. It is convenient to 
introduce the following notations: 



h = r]^"'h 
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and the we define the Christoffel symbols according to: 



The expression 

Rpv;pcj = dpT^,.y„ -{p^ a) (2.3) 

is called the Riemann tensor and it is gauge invariant for massless and massive gravity also. In 
the case of massive gravity we have new gauge invariants namely the (symmetric) tensor 

(ppu = -dpVi, - dyV^ + ripydpVp + m h^.^ (2.4) 

and its trace: 

</> = V""" 4>pu. (2.5) 
These expression are immediately proved to be gauge invariant. The same is true for their deriva- 
tives and the traceless part of these tensors. Let us denote by R^^l-p^-x^ \^,(k^pu-,pi...pn,(kfpi...pu 
the traceless parts of these tensors. We denote the co-cycles of dq by Zq. Then we have the 
following result [5]: 

Theorem 2.1 Let p S Zq. Then p is cohomologous to a polynomial in the traceless variables 
described above. 

We note that in the case of null mass the operator dq raises the canonical dimension by one unit 
and this fact is not true anymore in the massive case. We are lead to another cohomology group. 
Let us take as the space of co-chains the space P^"-) of polynomials of canonical dimension oj < n; 
then Z^Q^ C and B^q^ = dqT^""-^^ are the co-cycles and the co-boundaries respectively. It 
is possible that a polynomial is a co-boundary as an element of V but not as an element of "P*^""^ . 
The situation is described by the following generalization of the preceding theorem. 

(n) 

Theorem 2.2 Let p £ Z}^ ' . Then p is cohomologous to a polynomial of the form pi + dQP2 
where pi G Vq and p2 G 

We will call the co-cycles of the type pi (resp. dQP2) primary (resp. secondary). Using 
this result one can determine the most general form of the interaction between the massive 
gravity and a scalar field of mass M. We will call expressions of the type dqB^ + idpb^^ relative 
coboundaries. 

Theorem 2.3 Suppose that the interaction Lagrangean T between the massive gravity and a 
scalar field is trilinear in the fields (and their derivatives). Then T it is relatively cohomologous 
to the following expression: 



T = ci 



.^Ct)"" + C2$0^ + C3$^ </) + C4 (^dp^d^^h"" - ^M^<^^h^ + C5^>^ (2.6) 

i. e. 

dqT = i dpT'' (2.7) 

with 

= C4 (^^uf^d^^d^^ - u''d^'<^d^,^ - ^M^''^^^ . (2.8) 
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Proof: (i) By hypothesis we have ()2.7p and the descent procedure (based on a variant of the 
Poincare lemma [3]) leads to 

where the carets indicate antisymmetry and can choose the expressions 7^^^ to be Lorentz co- 
variant; we also have 

gh{T^) = \I\,uj{T^) <5. (2.10) 
^From the last relation in 12.91 we find, using the preceding Theorem 2.2, that 

J^ll^upa] ^ d^B^f^'^P^] + J^lf^'^P^] (2.11) 

with Tq^'^P'^^ G Vq^^ depending only on the invariants. It is easy to prove that such a (trilinear) 
expression does not exists so we have 

rplf^iypa] ^ dqE^'^P". (2.12) 
The third relation of the descent equations gives: 

dQ{T^P"P^ - i d^B^P'^P''^) = (2.13) 
so we obtain again with the preceding Theorem 

rplpup] ^ Q[pup] ^ . Q^Qit^Upcr] + r|.\^^^P] (2.14) 

where Tq''^''' £ depends only on the invariants. Again, we can see that such an expression 
does not exists so we have 

j^[pvp] ^ qIpM ^ • d^B^^^^^pcr]_ (2.15) 
The second descent equation then gives 

dQ{T^^"'^ - i dpB^^^'p'^) = 0. (2.16) 
(ii) We obtain from the relation (|2.16p with the preceding Theorem 2.2 

tI'^'^I = dg^I^^l + i dpB^^^'P^ + rd^'^l (2.17) 
where Tq'^'^^ G T^q^^ . The first descent equation gives the restriction: 

(iQ(r^ - apSt^^i) = a^TfS'"'^ (2.18) 

so the divergence diyT^P''^^ must be a coboundary. We do have a nontrivial expression for Tq^'^^ 
given by secondary cocycles. In the even sector with respect to parity we have 

tIp"^ = g^uPu''<^ + g2U^PP'^u^'"'^r]p^^ + gsu^P'^^uPdp'^ + giiu^PP^u" - u^''P'^uP)dp^ 

+g5{u^PP^d''<l> - u^''P^dP<l>)up + geiuPupd^dP - u^UpdPdP) (2.19) 
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and in the odd sector we have the expression e'^'^'"^?'[po-] where T^p^^^ has the same form as above 
but with Qj — > g'j. Here we have used the following notation: 

JH = Qi^^u _ Qv^i, ^2.20) 



One computes the divergence 5,t['^"1 and requires that it is a coboundary. After some compu- 
tations one finds out that the remaining terms can be grouped into a relative coboundary i.e. 
rj''''] = dob^"" - idph^^'^'py It follows that we have 

= dg^l'^'^l + i dpB^^"'P^ (2.21) 

if we redefine the expressions B^^'^^ and B^P'^f\ 
The first descent equation gives 

dQ{T^' -idpB^P''^) = {) (2.22) 
so if we use the Theorem 2.2 we find 

= dqB^" + iapB^^""^ + (2.23) 

where £ V^^\ If we substitute this in the starting relation (j2.7p we get the consistency 
condition 

dgiTP - tdpBP) = id^T^^ (2.24) 

i.e. the divergence O^Tq must be a coboundary. The generic form of Tq is again a secondary 
cocycle. In the even sector with respect to parity we have: 

= f^u^^^ + f2u>'d''<^du^ + fsu^d^^d^^ + Uu^'^'^d^^ + /5U^«>a^a^$. (2.25) 

In the odd sector we have 

= feP^'P'^u^p^d^^. (2.26) 
We compute the divergence dpT^ and the consistency condition leads to 

= f Qn'^r^S,^ - u''dP<^d,^ - ^M2n^$2^ + dgb^ + d.b^" (2.27) 

for some arbitrary constant /. One can get rid of the relative coboundary by redefining the 
expressions B^ and i?^*^. Moreover one proves that dpTQ = —idgt where 

t = f (hp^dP^d"^ - ^m2/i^>2 ) (2.28) 
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The starting relation (12. TP is now 

dqiT - t - idpB") = (2.29) 
so that a final use of the Theorem 2.2 gives 

T = t + dqB + idpBP + To (2.30) 

with Tq € Vq'^ . The generic form of Tq is 

To = Ci$4°V^°^'''' + C2^<A^ + C3$2 </> + C4$^ (2.31) 

The expression from the statement follows easily: we can replace by (pp,^ if we redefine the 
constant C2 and follows from (j2.27p . ■ 
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3 Second order gauge invariance 

In second order we must construct chronological products T{x,y) and T^{x,y) such that 

d 

dQT{x,y) = i—T^'{x,y) + X ^ y (3.1) 

is verified. The construction procedure is well-known: one first computes the causal commutators 
[T{x),T{y)] and [T^{x),T{y)] and substitutes the causal Pauli-Jordan distributions in the tree 
graph contributions by Feynman propagators D^{x — y). If on the right-hand side of (j3.ip a 
wave operator d"^ operates on we obtain a local term ~ 5{x — y). These anomalies must be 
compensated by finite renormalizations. 
The generic form of the anomaly is 

A{x, y) = 5{x - y)a{x) + [5^5(x - y)]a^'{x, y) (3.2) 

The total anomaly is obtained by adding the contribution A{y, x) with x, y interchanged. Then 
the terms with d5 can be combined by means of the identity 

[dl8{x - y)]/(x, y)+x^y=[dlf- d^^My - x)], (3.3) 

which follows by smearing with symmetric test functions; this is the right test function space 
here, due to the symmetry of the chronological products. Then the total anomaly is equal to 

Atotix, y) = [2a{x) + dla^' - d^X]5{x - y) = A{x)5{x - y). (3.4) 

The cancellation of the anomalies is equivalent to 

Atot{x, y) = dQR{x, y) - i9^i?^(x, y) + x ^ y; (3.5) 

here the expressions R{x,y) and R^{x,y) are finite renormalizations: these are quasilocal oper- 
ators: 

R{x, y) = 6{x - y)B{x) + ■■■ (3.6) 

and 

R^'{x,y) = 5{x-y)B''{x) + --- (3.7) 

where B and B^ are some Wick polynomials and • • • are similar terms with derivatives on the 
delta distribution. Indeed, in this case one can eliminate the anomaly by redefinition of the 
chronological products 

T{x,y) ^ T{x,y) + R{x,y) (3.8) 

and 

T^'{x,y)^T^'{x,y) + R^'{x,y). (3.9) 

One can prove that the cancellation ()3.5p of the anomalies is achieved if we can write the 
operator part ^(x) in (j3.4p in the form 

A = dQB-idf,B''. (3.10) 

In fact, the derivative terms in (13. 2p can be combined with help of the identity 

d^^iBf^ix, y)5{x -y)]+x^y = [d^^B'^ + dyBf^^x - y). (3.11) 
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The terms in which generate anomalies are the following ones: 

-A^u'^dav'^d^v^ - au^'da^df,^. (3.12) 

Here we have put the gravitational coupling constant k = 1 for simplicity. According to theorem 
the first order coupling to the scalar field <I> of mass M is given by 



+C4 ( d^^d^^h^"" - -M^^^h ] + C5^>^ (3.13) 



T$ = ci^cpf.u^b^" + C2$(m/i + 2df,v^'f + C3^'^{mh + 28, 

1 
2 

We first consider the couplings linear in $, i.e. with coefficients ci,C2. 

Theorem 3.1 Second order gauge invariance implies ci = and C2 = 0. 

Proof: To prove this result it is sufficient to find anomalies with ci or C2, which cannot be 
compensated. For ci we consider the commutator 

- Scu^dxv'id^v.ix), r^ivMapiy)^ 

(3.14) 

As described above the commutator gives a causal propagator which in the chronological product 
becomes a Feynman propagator. Applying the derivative d/dx'^ we get a d'^D^ leading to the 
anomaly 

Ai=Aiciu^dxv''{x)^y){2<l,Md^ -(t^{y)dl)5{x-y). (3.15) 

In the same way we consider the commutator 

- Sc2u''do,v''[d^v,{x),ct){ymy)^{y). (3.16) 

Here the resulting anomaly is equal to 

A2 = -^ic2u''d^v''{x)^y)<l,{y)dl5{x - y). (3.17) 

There are no other anomalies with Wick monomials uv<j)^, uv<j)^i,^, respectively. Consequently, 
Ai and A2 must cancel against each other. For the last Wick monomial we see from ()3.15p that 
ci must be and hence, C2 must also vanish. 

The situation is non-trivial for the remaining couplings which are bilinear in 

Theorem 3.2 Second order gauge invariance implies C4 = —2, but C3 and C5 remain unre- 
stricted. In the second-order chronological products the following finite renormalizations are 
necessary 

T{x,y)=T^{x,y) + i5{x-y)N{x) T^{x,y) = {x,y) + i6{x - y)N^'ix) (3.18) 

where 

N = 2$2/m2(2/i^^ V - h^) + C3 L(2/i^^ V - h^) + -{d^v^'d.v'' - ^^v^'^,v^') 
L L m 



-c^vi'd^^} (3.19) 
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and 

-2c3{2ui'dav'^ - u^dav")^'^. (3.20) 

Proof: 

In this proof we must calculate all anomalies containing We also give the commutators 
where the anomalies come from. From 

{-u'^do^h - dau'^h + 2d''u''K,)[d^h{x), h{y)] imc3 - —cA 

we get the anomaly 

Ai = 2i{2mc3 - M^C4){-u''dah - dau'^h + 2d''u''hau)^^S, (3.21) 

and 



leads to 

A2 = i{2mc3 - M'^ci){u^dxh + d^u^h - 2dxu"'Kx)^'^5, (3.22) 

The commutator 

{-u'^dah - dau'^h + 2d''u'^hau)[di^h{x), hf^^{y)]c4dp^d^^ 

gives 

A3 = icA{-u'^dah - dau'^h + 2d''u'^Kv)dp^dl^^ (3.23) 

and 

2{u^dxh'"' + dxu^h'^'' - dxu'^h^'' - dxv^h'^^)[d^Ky{x), /i^^(y)]c45/3$5^$ 

yields 

A4 = —iC4 

-2dxu^h?^ - 2dxu'h^^ + 2dxuah^''r]'^^]c4d(3^d^^. (3.24) 
Next the commutator 

- 4u''dc,v''[di,v,{x), dpvf^{y)]2c3^'' 

leads to 

As = -4ic3u''dav''{x)^''{y)dyS{x - y). (3.25) 

and finally 

- C4U^dp^ [d^^^ix), $'(y) {mc3h + 2c3do,v'^ - ^M^h)j + C4daHy)d,^h'''' + C5^^{y) 



gives 

Aq = iciuf^dp^{x){2^ (mc^ih - yAf^) h + Acs^d^v" + 2cida^{y)h'"' {y)dl + 3c5$^}(5(x - y). 

(3.26) 

The sum Ai + . . . + \s equal to 

Bi = -2ici{u^dxh''^ + ^AU^/i"^ - ^AM^/i"^ - dxu''h^^)da^dp^5 

+2iclu'^d(i^{x)h^'''{y)d^^{y)dl5{x - y) (Tl) 
-i(2mc3 - M'^C4){uf'd^h + 

+ic4(2mc3 - m2c4)u'^5/3$/i$ (r2) 
+2z(2mc3 - M^C4)d^u^h'"'^'^6 (TS) 

+3ic4C5M^9/3$$^(5. (r5) (3.27) 

Following the methods developed in [9j (Sect. 5. 9) we have grouped the terms according to their 
type of Lorentz contractions. For example, (Tl) has u^h°'^^^ and 3 derivatives which is different 
from (T3). Only the terms within one type Tl, . . . T4 can be combined to give a divergence. Due 
to the different coefficients C4 and c| in Tl we must have C4 = —2 in order to get a divergence. 
If C4 were 7^ —2 then the last term of (Tl) would remain without compensation. Since this term 
is not a relative coboundary gauge invariance then would be violated. 

The total anomaly is obtained by adding the contribution x <^ y according to ()3.ip . For the 
terms with 5{x — y) this simply gives factor 2. For the terms with derivative of 5 we use the 
identity 

g{x)f{y)dl5{x - y) + x ^ y = {8^9 f - gdaf)S{x - y) (3.28) 
Now the total anomalies of type Tl in (3.14) can be written in the form 

(ri)tot = -4ic4 [{u^dxh^^ + dxu^h"^)d^^dp^ - dxu^h^^d^^dp^ 

+u'^dad,3^dx^h''^ - uf^d/B'^da^dxh"^ - df3<^dadx<^>h''^ . (3.29) 
This agrees with the result in massless gravity [9], eq. (5.9.40), and is a divergence 

(ri)tot = -^icAdf [(n^/i°^ - u'^h^'^)da^dp<l>6{x - y)] 

+x ^ y, (3.30) 

where C4 = —2 has been taken into account and will be assumed in the following. Type T2 is a 
divergence as well: 

(T2)tot = -2i(M2 + mc3)a^K/i$2^(x -y)]+x^y. (3.31) 

As in the massless case ([9], eq.(5.9.45)) T3 is a coboundary: 

(r3)tot = 2(m2 + mc3)dQ[{h^ - 2h^,h^"')^^5{x - y)]. (3.32) 

Using the identity (j3.28p we write T4 as follows 

(T4)tot = -4ic3[(5^u^ + u^'d^d^v'')^^ - 2u^'df,v''<^d^,<^ 

+Au^'d^v''^di,^]6{x -y). (3.33) 
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We first split off a divergence 



(r4)tot = -4ic3[29^(n'^9,7;^$2) _ d,{u''d^v''<^^) 

^1{d^u^d^v^ - ^^,u^^^v'')^'^\b{x - y). (3.34) 

Now tlie terms in tlie second line are a coboundary 

(T4)tot = -2ic3[2d^^{u^'d,v^<l>H) - d:u''d^v^^^5)]+x^ y 

+8^dQ[(5,^;^9X - d^vf'd.v'')'^^. (3.35) 



Finally, T5 is a coboundary 

T5 = —C5dQ{vf'<^^df,^S). (3.36) 
m 

Adding the contribution x <-> y this gives the result of the theorem. ■ 

4 Modified general relativity 

As in ref. [1] we now consider the limit m ^ of vanishing graviton mass. The point is that 
this does not lead to massless gravity because the vector graviton field v^^ does not decouple 
from the other fields. In fact, in first order (proportional to Newton's constant) there survives 
the coupling term 

= h^^'d^d^vx. (4.1) 

If scalar matter is included then in addition to the standard coupling to the energy-momentum 
tensor of the scalar field (~ C4 in (13.13P ) two further couplings T3 and T5 are possible. However, 
in second order the graviton mass appears in the denominator in N in (j3.19p . Consequently, if 
the coupling constants C3 and C5 do not depend on m, the limit m ^ exists for C3 = = C5, 
only. Then we have no direct coupling of the ii-field to normal matter; the resulting theory of 
[1] seems not to be physically relevant. 

There is another option. Gauge invariance does not forbid the possibility that C3 and C5 
depend on m, for example 

C3 = Asm, C5 = Asm, (4.2) 

where \j are independent of m. Then in the limit m — > the first order trilinear couplings die 
away, but there remain the following quartic couplings from second order 

= 2$2|8A3(a^t;'^9.t;" - ^,,v''^,v^') - 12A5t;^a^«>}. (4.3) 

In the classical limit this coupling must be added to the classical Lagrangean. As in the other 
coupling terms the usual factor \/—g is included. The necessity of this factor becomes clear 
when we derive the field equations below; but of course, an independent check by a third order 
calculation must be done. Our modification of general relativity is now defined by the following 
Lagrangean density 

Ltot = ^V^R + V^g'^'d^vxd^v^ 

+lv^g'"'d^^d,'^ + V^^^[x3{d^v''d,v'' - d^v^d^v^ + X5vf'd^i>]. (4.4) 
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The two terms in the first line are the pure gravitational interactions which have been studied 



already in [T]. The first term is the Einstein-Hilbert Lagrangean, 



.-2 



327rG is essentially 



Newton's constant and R the scalar curvature. The second line contains the interaction with 
scalar matter; the numerical factors in (j4.3p have been absorbed by redefining the coupling 
constants A3 and A5. 

The Lagrangean ()4.4p as it stands is Lorentz invariant, but the new terms in the second 
line are not invariant under general coordinate transformations. In [l] we have argued that this 
latter invariance can be maintained in the second term of the first line, if we consider as 
four scalar fields. This argument cannot be used for the new matter couplings in the second 
line. The lack of general covariance might be disturbing for classical relativists. However, one 
should keep in mind that classical general covariance corresponds to gauge invariance of the 
spin-2 quantum gauge theory, so that this principle is incorporated in the quantum theory. The 
latter is background dependent; we have selected Minkowski background. Returning again to 
the classical theory this background dependence remains and we get a Lorentz invariant classical 
theory, not a general covariant one. Still, by checking gauge invariance in third order we have to 
test whether there are further modifications in the classical theory. This will be done elsewhere. 

The Euler-Lagrange equations for the Lagrangean (j4.4p give the system of coupled field 
equations. Variation of g'^'^ gives the modified Einstein equations 



R 



The variational derivative with respect to yields 

dai^gg'^^dpy^) = 2A3 [d^i^g^^d.v^) - d.i^g^^d^v'') 



(4.5) 



(4.6) 



Here the vector-graviton field has source terms from the new scalar-matter coupling. Note 
that the second order derivative dfj,di,v^ cancels on the right-hand side so that we have a wave 
equation with source. Finally, the variation of $ gives the Klein-Gordon equation in the metric 
gO!-P p^^g source terms: 



1 



--d„ 



-9 



(4.7) 



The physical consequences of these field equations remain to be investigated, in particular 
whether there are solutions giving an explanation of the dark matter phenomenology. 
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